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ABSTRACT 

The feedback from galactic outflows is thought to play an important role in shaping the gas content, 
star formation history, and ultimately the stellar mass function of galaxies. Here we present evidence 
for massive molecular outflows associated with ultra-luminous infrared galaxies (ULIRGs) in the 
coadded Redshift Search Receiver 12 CO (1-0) spectrum. Our stacked spectrum of 27 ULIRGs at 
z = 0.043 - 0.11 (z/rest = 110 - 120 GHz) shows broad wings around the CO line with AV(FWZI) w 
2000 km s _1 . Its integrated line flux accounts for up to 25±5% of the total CO line luminosity. When 
interpreted as a massive molecular outflow wind, the associated mechanical energy can be explained 
by a concentrated starburst with SFR > 100 M Q yr _1 , which agrees well with their SFR derived 
from the FIR luminosity. Using the high signal-to-noise stacked composite spectrum, we also probe 
13 CO and 12 CN emission in the sample and discuss how the chemical abundance of molecular gas may 
vary depending on the physical conditions of the nuclear region. 

Subject headings: galaxies: evolution — galaxies: interactions — galaxies: ISM — galaxies: nuclei — 
galaxies: starburst — ISM: jets and outflows 



1. INTRODUCTION 

Extreme star formation rates of > 10 3 M® yr" 1 
have been derived for lumin ous infrared galaxies dis- 
covered by deep IR s urveys (|Downes fc Solomonl 119981 : 
lYounger et all l2008t iRiechers et al.l 12009ft . In such 
radiation-pressure s upported galactic disk s with the 
maximum starburst (|Thompson et al.ll2005D . large scale 
outflows can be triggered by superwinds from mas- 
sive young stars and supernova explosions, which may 
play important role s in galaxy formation and evolution 
ICooper et al.l (|2008f ) . Galactic outflows can regul ate star 
formation by heating cool gas (|Tang et al.l l2009h . They 
can enrich both interg alactic medium and galactic disks 
(|Heckman et al.lll990l) . Their feedback can also explain 
the apparent discrepancy between the theoretical predic- 
tion of the dark matter halo mass function and the mea- 
sured stellar mas s function for galaxies in th e successful 
ACDM scenario (|Springel fc Hernauistl[2003l ). 

In this Letter, we probe kinematic signatures of molec- 
ular outflows in a sample of 27 ultra-lumi nous infrared 
galaxi es (ULIRGs) recently studied by IChung et al.l 
(2009). A large fraction of UL IRGs are mainly pow ered 
by merger- induced starbursts ([Sanders et al.lll988T ) and 
hence make good targets for investigating associated out- 
flows. Evidence for such an outflow has been reported in 
a sing le ULIRG syste m such as Arp 2 20 (Saka moto et al.l 
120091 ) and Mrk 231 (jFeruglio et all IMlbT) . Those out- 
flow signatures are however too faint to be detected in- 
dividually in our ULIRG sample. Therefore we employ a 
stacking analysis to look for faint and broad high velocity 
line wings in the 12 CO line profile. In order to inspect 
outflow driving mechanisms, our ULIRG sample is par- 
titioned into two groups based on optical emission line 
diagnostics: starburst dominated galaxies and galaxies 
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with large AGN contributions. With the reduced noise 
in the stacked composite spectrum, we also measure the 
average brightness of other weaker molecular lines such 
as 13 CO(1-0) and 12 CN(l-0), as an independent mea- 
surement of molecular gas properties, that have been 
detected only in the ne arest IR luminous galaxies (e.g. 
lAalto et al.iri995ll2002h . 

2. SAMPLE AND STACKING 

We use the sample and the data from the recent Red- 
shift Search Receiv er (RSR) 12 CO J = 1 — > survey 
of local ULIRGs by IChung et all (|2009ft . The observa- 
tions were carried out with the Five College Radio As- 
tronomy Observatory (FCRAO) 14m Telescope in 2007 
and 2008, targeting 29 ULIRGs at z = .043 - 0.11. As 
discussed in detail by IChung et al.l (|2009h . this is a repre- 
sentative subset of ULIRGs as the primary selection cri- 
teria were those related to observational scheduling and 
the redshift range that brings the 12 CO J = 1 — > line 
within the bandpass of the RSR system. In our stack- 
ing analysis, we include only the 27 CO detected objects. 
The CO line luminosity L' co of the sample ranges from 
1.2 to 15.3xl0 9 K km s _1 pc 2 with a median value of 
6.7xl0 9 K km s" 1 pc 2 . 

The stacking of the RSR spectra is performed using 
the following procedure. First, each coadded spectrum is 
shifted to the rest frequency by multiplying the observed 
frequency by (l+z co ), where z co is th e CO redshift of eac h 
ULIRG derived from the line fitting ()Chung et al.ll2009ft . 
Each CO spectrum is normalized by the best fit Gaussian 
peak, and then all spectra are "aligned" at the frequency 
centroid. A linear interpolation is used in the alignment 
process. The normalized spectra are averaged, weighted 
by the rms "noise" measured in the normalized spectra, 
excluding the ±0.5 GHz regions around the three tran- 
sitions of our interest, 13 CO, 12 CN, and 12 CO as well 
as the noisy end channels. Finally, the averaged spec- 
trum is Hanning smoothed to produce the final spectral 
resolution of 61 MHz (158 km s" 1 at 115.27 GHz). 

A "non-ULIRG" comparison spectrum was derived by 
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stacking the RSR spectra of 19 z = 0.037 — 0.066 galax- 
ies selected for their high HI mass (Mm >2x 10 10 M ; 
Haynes et al. in prep., O'Neil et al., in prep). These 19 
galaxies were selected from another RSR commission- 
ing programs - the CO survey of 29 HI rich galaxies 
at similar redshifts. Nineteen galaxies were detected in 
CO with comparable S/N and sensitivity as our ULIRG 
sample, but they are otherwise normal in their star for- 
mation and nuclear activities. These HI rich galaxies 
mostly look like normal spirals in the optical and their 
mean FIR luminosity is 2.8 ± 1.4 x 10 10 L Q , 30 times 
lower than that of our ULIRG sample. Their L' co of 0.4- 
3.2 xlO 9 K km s _1 pc 2 is only slightly smaller than that 
of our ULIRG sample. Further details of the RSR CO ob- 
servations of these Hi-rich spirals will be presented else- 
where (Chung et al. in prep). 

The rms noise in the final spectra for the ULIRG sam- 
ple and the control sample, normalized by the 12 CO peak 
flux, are 0.014 and 0.027, respectively, yielding the S/N 
that is better than those of individual spectra by a fac- 
tor of 5 to 47 (see Table [TJ. Figure [1] shows the stacked, 
normalized spectrum of a broad frequency range (109.65- 
116.25 GHz) which includes all three transitions, 13 CO 
(1-0), 12 CN (1-0), and 12 CO (1-0). In Figured we zoom 
in the 4500 km s" 1 range around the 12 CO line to show 
the characteristics of the profile more in detail. 

3. STARBURST POWERED OUTFLOWS 

The stacked spectrum of the ULIRG group has re- 
vealed broad wings around the CO line as seen on the 
top of Figure [TJ. The wings are blue- and redshifted 
by 1000 km s" 1 from the main CO line peak (FWZI 
of 2000 km s" 1 ) with the total line integral of 19±5% 
of the total (see Fig. |2J. The line wings of the ULIRG 
stacked spectrum are detected with S/N~ 3 in each chan- 
nel, which would be difficult to be detected in individual 
spectra. In fact, the effective integration time of the 
stacked ULIRG spectrum is 115.8 hrs, 10-60 times more 
than the integration on individual ULIRGs. Note that 
such wings are not present in the control sample of HI 
rich galaxies with the effective integration 127.7 hrs as 
shown on the bottom of Figure [TJ The comparison of 
the ULIRGs with the non-ULIRG population is better 
shown in the upper two panels of Figure [5] 

Such broad wings can form when entrained cool 
gas gets ejecte d alo ng with hot ioniz ed outflowing gas 
(|Curran et all 119991; iNaravanan et al.ll2006D . In order to 
examine whether a starburst or an AGN is powering the 
outflow, we have divided the ULIRG sample into two 
groups: (1) 14 ULIRGs which are classified as "Sbrst" or 
"Hn" (SB group) with no obvious sign of AGNs; and 
(2) 13 ULIRGs with "seyfert" spectra (AGN group). 
This grouping is done based on the classification from 
the NASA/IPAC Extragalactic Database (NED) 4 . Ob- 
jects with a "LINER" classification are included in the 
AGN group unless it also has a "Sbrst" or "Hn" des- 
ignation. Galaxies with a hybrid ("Sy+SB/Hn") clas- 
sification are included in the AGN group. Among our 
sample, eight ULIRGs from each group ha ve been mod- 
eled in their spectral energy distribution bv lFarrah et al.l 

4 This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 
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Fig. 1. — Full Composite Redshift Search Receiver (RSR) spectra 
of the ULIRG and comparison sample. The rest frequencies of 
13 CO (110.20 GHz), 12 CN (113.50 GHz), and 12 CO (115.27 GHz) 
are indicated with dotted vertical lines. The stacked spectrum 
of 27 12 CO detected ULIRGs is shown on the top, the stacked 
spectrum of a subsample of 14 "Sbrst" or "Hn" (SB group, i/^f N < 
0.12L IR ) and 13 seyfert or LINER type ULIRGs (AGN group, 
L*^ N ss 0.32Ljji) are shown in the middle, followed by the stacked 
spectrum of 19 non-ULIRG sample at the bottom. The FWZI 
spectral regions used to measure the line and wing flux density are 
shown as boxes. 



(2003) who found an AGN contribution to the IR lumi- 
nosity of >27% for the eight in the AGN group, which 
is a factor of two higher than that of the eight in the 
SB group (~ 12%). The mean S/N measured by the ra- 
tio of the CO line peak and the rms is highest for the 
SB group due to a few objects with the strongest CO 
emission among the sample, but this ratio does not vary 
significantly from group to group. In fact, the mean rms 
of each group and the CO linewidths of different groups 
are very similar as summarized in Table [TJ making our 
results robust. 

We show the comparison of the two groups on the bot- 
tom of Figure^] The rms noise in the normalized stacked 
spectra is 0.0189 and 0.032 for the SB group and the 
AGN group, respectively. The wings around 12 CO line 
appears to be even stronger in the SB group compared to 
the wings seen in the entire ULIRG sample, with more 
flux in the wings which is about 25% of the total CO flux. 
These broad features however disappear when only the 
ULIRGs with Seyfert spectra are combined. The total 
line flux was measured by integrating the line flux density 
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TABLE 1 

Outflow Properties and Molecular Abundance 



Group 



Moan 

rms (mK/TX) T peak /rms W co (km/s) 



Stacked 
12CO 



ULIRGs 

All (27) 0.072±0.023 0.46±0.15 6.4±4.1 263± 59 0.19±0.05 

Sbrst+Hn (14) 0.071±0.027 0.45±0.15 7.3±4.2 266± 67 0.25±0.06 

Sy+LIN (13) 0.073±0.026 0.47±0.16 5.4±4.0 262± 54 < 0.12 

Non-ULIRGs 

Hi-rich spirals (19) 0.050±0.007 0.39±0.09 5.0±2.7 266±103 < 0.11 > 7.3 > 7.3 



> 16.6 > 16.6 

> 13.3 > 13.3 
11.1±6.2 9.3 ± 4.3 
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Fig. 2. — A zoom-in view of the composite CO J = 1 — > 
Redshift Search Receiver (RSR) spect ra. On the top, the ULIRG 
sample of 27 CO detected galaxies from Chung ct al. (2009) is com- 
pared to the non ULIRG sample of 19 Hi-rich galaxies. On the 
bottom, the CO spectra of the SB group and the AGN group are 
compared. 



within the FWZI regions (thin solid lines in Fig. 1) for 
both groups with and without wings. The wing-only flux 
was measured within the same FWZI as the AGN group 
and the control sample, and the difference between the 
total and the line wing flux has been adopted as the CO 
flux for these groups. The fractional flux in the wings and 
the upp er limits of different subgroups are summarized in 
TableQ] iRupke et~a l. (2005) also found a lower frequency 
of neutral wind among Seyfert 2 ULIRGs in their study 
of Na I D absorption line in 26 AGN/starburst-composite 
ULIRGs at 0.03 < z < 0.44, further supporting the star- 
burst origin for this massive neutral wind. 

The energetics of the neutral wind traced in broad CO 
wings is also consistent with being powered by the on- 
going starburst traced in the far-infrared. The energy 
injection rate (dE/dt) in a wind-blown bubble of a ra- 
dius R expanding at velocity v into an infinite homo- 
geneous medium w ith density tlq can be expressed as 
(Weaver et allll977ft . 

dE 

— - 3.3 x 10 35 Rl pc v 3 km/s n . cm -3 erg s" 1 . (1) 
Adopting a bubble size of 0.2 kpc lSakamoto et al.l ([2009D 



found for the high velocity molecular wind in Arp 220, 
an outflow velocity of 1000 km s _1 from the line wing 
velocity in the stacke d spectrum, and an ambient den- 
sity of 10 cm" 3 (e.g. IVeilleux et al.lll995h . we drive an 
energy injection rate of ~ 1.3 x 10 44 ergs s _1 . As- 
suming an energy ou tput per supernova of ~ 10 51 ergs 
(jVeilleux et al.lll995l ). our estimated dE/dt yields a su- 
pernova rate, ^sN.yr- 1 °f 4 yr . Using a Scalo initial 
mass function (IMF) with a mass range of 5 — 100 Mq, 
the star formation rate inferred from this supernova rate 
(5F^m>5Mo = 24.4 ^sN,yr-i Mq vr" 1 . iCondonl [1991: 
lRosa-Gonzale3 12001 is SFR « 100 Mq yr" 1 . This 
agrees well with the SFR derived from the FIR lumi- 
nosity for these ULI RGs, 134-352 M yr" 1 (see Eq. 8 of 
iHopkins et aDl2003f ). 

The outflow speed implied by the CO line wings, 
1000 km s" 1 , is comparable to the wind velocity mea- 
sured by other phases of the superw i nd su ch as Ha 
and Na D f4 00-800 km s" 1 . iMartinl Il999l ) and OH 
(1400 km s~ 1 . lFischer et al.ll2010ft If thev all trace the 
same outflow, then Eq. (1) suggests that the spatial ex- 
tent of Na D and Ha should be larger (~ 0.4 kpc) while 
the spatial scale for the OH winds measured by Herschel 
would be more compact (< 0.1 kpc) than the CO wind. 
The fact that wind velocities measured in molecular out- 
flows are larger than in the optical may indicate that 
supernova driven wind embedded in molecular gas slows 
down by the time it breaks out of the starburst region. 

The inferred outflow speed exceeds the escape velocity 
and is high enough to blow away the molecular gas that 
hosts the star formation activity and to pollute the sur- 
rounding IGM significantly. Depending on whether the 
CO line is optically thin or thick, the outflowing molecu- 
lar gas mass ranges between 1-6 x 10 9 Mq. The line wings 
are symmetric in intensity and shape on both sides of the 
line, and this suggests that the wings are bipolar in geom- 
etry and likely optically thin. This is not contradictory 
to the highly asymmetric CO 3 — > 2 line wings in Arp 220 
(FWZI = 1000 km s" 1 ) found bv lSakamoto eTall (|2009n 
since the J = 3 — > 2 line has a higher optical depth. 
These observations suggest that more than 10 9 Mq of 
molecular gas can be removed from the central starburst 
region through such a wind, rapidly depleting the gas 
supply for the starburst. Some of this gas may eventu- 
ally rain back on to the galaxy, enriching the galactic disk 
(|Heckmanll2003T) . 

Our conclusion that the central starburst can power 
th e observed massive o utflow contradicts the conclusion 
bv lFeruglio et al.1 (|2010f ) that the 750 km s^ 1 wind in CO 
1 — > found in Mrk 231 is powered by the AGN activity 
and thus is an example of a "quasar feedback" at work. 
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Although iFeruglio et al.l (|2010t ) adopted a smaller than 
Galactic CO-to-H2 conversion factor, they may still have 
over-estimated gas mass and mass outflow rate as an op- 
tically thin estimate leads to a > 2 times smaller gas 
mass, lowering the mass outflow rate much closer to the 
current SFR. A stronger case for an AGN-driven molec- 
ular outflow is found in NGC 1266 where the optically 
thin CO mass outflow rate clearly exceeds the observed 
current SFR (Alatalo et al., in prep). The CO outflow 
velocity is much lower (~ 400 km s _1 ), however, and 
this phenomenon may not be very common, at least at 
low-z, since this is the only object with an AGN-driven 
molecular outflow found i n their survey of a large num- 
ber of early- type galaxies. iNaravanan et al.l ()2008f ) have 
shown that an AGN-driven molecular outflow may per- 
sist longer than a SB-driven outflow using a numerical 
simulation, but the accuracy of such model predictions 
and the sub-grid physics included needs to be tested fur- 
ther using a large sample of AGN+SB systems. 

4. 13CO AND 12CN ABUNDANCES 

Two important molecular transitions also appear in 
our stacked composite spectra, and we examine their 
line strengths to gain further insights into the molecu- 
lar ISM in these ULIRGs. Those are the lowest transi- 
tions of 13 CO and 12 CN (CN hereafter) at 110.20 and 
113.50 GHz, which have been detected in local starburst 
and Seyfert galaxies (ICasoli et all Il992l ; lAalto et"aH 
119951 l2002UPerez-Beaupuits et al.ll2007ft . 

In the RSR composite spectra, we find > 3 a bumps at 
both 13 CO and CN frequencies only in the AGN ULIRG 
group with the flux ratios, 11.1 ± 6.2 and 9.3 ± 4.3 for 
12 CO/ 13 CO and 12 CO/CN, respectively (Fig. [J). The 
other groups do not show such features, and the lower 
limits in 12 CO/ 13 CO and 12 CO/CN are summarized 
in Table [TJ The same linewidths as the AGN group 
have been adopted to calculate the upper limits of 13 CO 
(660 km s _1 ) and CN (600 km s _1 ) for the other groups. 

The ratio 12 CO/ 13 CO (R w ) has been reported to 
be generally larger i n starburst galaxies (Rio > 20, 
iGlenn fc Hunte r 2001) than in optically thick normal spi- 
rals (10 < R w < 20. ICasoli et al.lll992f) or Sevfert galax- 
ies (i?io ~ 12. iPapadopoulos fc Seaqu"istlll998l ). It has 
been suggested that the ove rproduction of 12 C, a p rimary 
product of nucleosynthesis (|Boreiko fc Betd ll996) in ac- 
tively star forming galaxies, is respo nsible for thi s trend 
(ICasoli et al.l Il992h . Alternatively, lAalto et al.l (|199lL 
also 1995) have suggested that i?io, which gauges the op- 
tical depth of * 3 CO gas in LTE (Ii2 CO /Ii3 CO ~ l/ Tl3 coi 
iPaglione et aLll2001h . can increase when molecular clouds 
are disturbed by powerful tidal force in merger driven 
starburst galaxies. Increased velocity dispersion within 
GMCs and a broader cloud-to-cloud velocity distribution 
can reduce the 13 CO opacity within these starburst nu- 
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clei. 

Meanwhile, CN is known to be a tracer of dense gas 
wi th lower critical density than HCN (by a fac tor of 
5; iPerez-Beaupuits et al.l l2007t iBaan et al.l 120081) . CN 
molecule is a p hoto- or X-ray d issociation product of 
HCN and HNC (jBaan et al.ll2008| ). a nd is predicted to b e 
abundant in both PDR s and XDRs (|Kohno et al.ll2008l) . 
iMeiierink et al.l (|2007l ) however found the CN/HCN ra- 
tio to be enhanced in XDR than in PDR, and toward the 
DR edges where the gas is highly ionized as in XDR. Our 
finding of the lowest CO/CN ratio in the AGN ULIRG 
group may imply a higher ionization rate of dense molec- 
ular gas as predicted by the Meijerink et al. models. 

We have only 8 (four SB and four AGN) and 14 ob- 
jects (nine SB and five AGN) whose 13 CO and CN lines 
fall within the RSR frequency band, and the significance 
of these results will have to be confirmed with a larger 
sample. 

5. FUTURE PROSPECTS 

There are ongoing theoretical efforts to model galac- 
tic outflows to understand their detai led properties such 
as th eir frequency and energetics (e.g. IChoi fc Nagaminel 
|2010[) and the feedbac k on scaling relations of galaxies 
(e.g. ISales et al.ll2010t ). Even for objects at cosmologi- 
cal distances where more direct morphological clues such 
as superbubbles, filaments, and chimneys are not visible, 
these outflow models can be tested by examining their 
spectroscopic signatures. Presently there is no consensus 
to on the drivin g mechanism for t he ob s erved outflows: 
iSakamoto et al.l ( i.e. s tarburst - I2009D ; iRiechers et al.l 
(i.e. starburst - I2009L vs. AGN - Ferglio et al. 2010; 
Fischer et al. 2010). Obtaining a better understanding 
is a pre-requisite in evaluating the importance of outflow 
feedback plays in galaxy evolution. The Redshift Search 
Receiver on the Large Millimeter Telescope (LMT) - a 
50m single-dish facility being built at Volcan Sierra Ne- 
gra, near Puebla, Mexico will extend our capability to 
study galaxy outflows and winds at higher redshifts with 
its vastly improved sensitivity. Spatially resolved mor- 
phological and kincmatical details obtainable using the 
ALMA will offer us the most stringent observational test 
for the origin of these massive molecular outflows. 
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